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ABSTRACT 



We present the discovery of a Baldwin effect in 8 nearby Seyfert galaxies for the three most prominent mid-infrared forbidden emission 
lines observable from the ground that are commonly found in AGN, [Arm] (A 8. 99 yum), [S rv](A 10.51 /jm), and [Nen](/112.81 /jm). The 
observations were carried out using the VLT/VISIR imager and spectroraph at the ESO/Paranal observatory. The bulk of the observed 
line emission comes from the inner <074 which corresponds to spatial scales <100pc in our object sample. The correlation index 
is approximately -0.6 without significant difference among the lines. This is the strongest anti-correlation between line equivalent 
width and continuum luminosity found so far. In the case of Circinus, we show that despite the use of mid-infrared lines, obscuration 
by either the host galaxy or the circumnuclear dust torus might affect the equivalent widths. Given the small observed spatial scales 
from which most of the line emission emanates, it is unclear how these observations fit into the favored "disappearing NLR" scenario 
for the narrow-line Baldwin effect. 
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1. Introduction 

The Baldwin effect, a negative correlation between equivalent 
width, Wa, and continuum luminosity, L, is commonly found 
in broad UV/optical emission li nes of AGN. In itially, it was re- 
ported for the Civ(^1549) line dBaldwinlll977l) for which W A oc 
■ L -0 - 2 , but later also found in emission lines of many other species 
for which the strength of the correlation d epends on the ioniza - 
tion potential of the respective line (e.g. iDietrich et al1l2002l) . 
Although the Baldwin effect is well established, its physical ori- 
gin is still a matter of debate. A common explanation refers to the 
spectral shape of the ionizing continuum which appears softer in 
high luminosity sources than in AGN with lower luminosities. 
This results in a relatively lower number of ionizing photons in 
the broad-line region (BLR) which, in turn, leads to a smaller 
Wa- While this explanation is, at least, qualitatively in agreement 
with the observed dependence on line ionization energy, it is not 
clear if L is the fundamental parameter that drives the broad- 
line Baldwin effect (BLBE). Instead, a similar effect has been 
observed if L is replaced by the bl ack hole mass , Mbh, or the 
Eddington ratio, l E dd = i/^Edd (e.g. lWande3ll999tlWarner et alj 

Hjoi. 

In addition to the BLBE, several narrow optical emission 
lines have been found to show a Baldwin effect with simi- 
lar order-of-magni tude luminosity-dependence as well (here- 
after: oNLBE: e.g.lBrorsen & GreerJ1992tlMcIntosh et aljl999t 
ICroom eFa l. 2002). Given the different spatial scales where the 
broad- and the narrow-lines are emitted, it is generally believed 
that the driving physical mechanisms for the BLBE and oNLBE 
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are different (e.g. IShieldsl 120071 and references therein). A 
luminosity-dependence of the equivalent width, sometimes re- 
ferred to as Iwasawa-Taniguchi effect, is also known for the 
narrow Ka-line in X-rays (e.g. Ilwasawa & Taniguchi 1993; 
iBianchi et al1l2007l) . 

In this letter, we present our discovery of a strong Baldwin 
effect in nearby Seyfert galaxies for the most prominent nar- 
row emission lines that are present in AGN mid-infrared spectra 
and observable from the ground, in particular [Arm](/18.99 yum), 
[Srv](A10.51//m), and [Nen](A12.81 jum). The Baldwin Effect 
in these mid-infrared forbidden lines (hereafter: iNLBE) is much 
stronger than for any other as yet observed. In Sect |2l we re- 
port on details about the VLT/VISIR observations. In Sect. [3] 
the results are presented, followed by a discussion of our find- 
ing in Sect. [4] This paper (hereafter: paper I) presents first re- 
sults on line emission of our larger campaign using VLT/VISIR 
which is dedicated to mid-infrared spectroscopy of nearby AGN 
at high spatial resolution. A forthcoming paper (hereafter: pa- 
per II) will deal with the mid-infrared continuum emission of 
AGN. In the following, we will use cosmologic al parameters 
H = 73km/(sMpc), Q A = 0.72, and Q„, = 0.24 dSpergel et al.l 
2007). 



2. Observations 

We used the mid-infrared imager and spectrograph VISIR lo- 
cated at the 8 m UT3 telescope of the ESO/Paranal observatory. 
With VISIR, we observed a sample of 8 nearby AGN of which 
5 are Seyfert 2 (Slh, Sli, S2) and 3 are Seyfert 1 (S1.0, SI. 5) 
galaxies. The 8 - 13yum A^-band spectra were acquired in low 
spectral resolution mode (R~300), for which 4 different wave- 
length settings centered at 8.5, 9.8, 11.4, and 12.4//m have to 
be combined for each object. With VISIR, the achieved spatial 
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Fig. 1. Equivalent width, Wa., of the [Arm] (blue circles) and 
[S iv] (red crosses) lines plotted over the X-ray luminosity, Lx, 
for our sample of nearby AGN as listed in Table Q] The fit- 
ted correlations Wa k L x 075±011 for [Arm] (blue dashed line; 
pspearman = -0.9 1, significance 0.002), and W A oc L^>.66±0.12 for 
[S iv] (red dashed-dotted line; pspearman = -0.79, significance 
0.04) are shown. Circinus (gray symbols) is outlying probably 
due to significant obscuration by dust in the host galaxy or the 
circumnuclear dust torus (see Fig.|4]i. 



resolution in this wavelength range is 0725 - 0739, which is a 
factor 2-4 smaller than the used slit widths of 0775 - 170 mak- 
ing slit losses negligible. The resulting spectra were extracted 
using the standard VISIR CPL pipeline (V3.0.0) and calibrated 
by a number of generic standard stars. For NGC 4507, a single 
standard star observed within the same night has been used for 
calibration. In addition to the 8 objects, archival VISIR data of 
NGC 1068 and NGC 253 has been downloaded, extracted, and 
calibrated in the same way. A more detailed description of the 
observations will be presented in paper II. NGC 4593 has been 
observed only in the 8.5 and 9.8 //m settings, corresponding to a 
restframe wavelength coverage of approximately 7.9 - 10.4 //m. 
As a result, only [Ar m](/I8.99 ^m) is seen in the combined spec- 
trum, giving Wa data in this line on a total of 8 objects. For 
[S iv](/I10.51 jum), and [Nen](/I12.81 //m), equivalent widths of 
7 objects is available. 

3. The Baldwin effect in forbidden mid-IR lines 

3.1. Results 

In Table[U we present extracted Wa of all three lines. For all ob- 
jects, the line emission was strongly peaked towards the nucleus. 
Except for the emission in the most nearby objects NGC 1068 
(~ 10% extension in [Siv] with respect to the continuum PSF) 
and Circinus (~50% extension in [Siv] and [Nen]), no spatial 
extension was detected in the continuum and the lines. The re- 
duction was done in a way that even for these extended objects, 
almost no line flux is lost by the extraction window. The errors 
on Wa represent the uncertainty due to the determination of the 



Fig. 2. Equivalent width, Wa, of the [Nen] line plotted over the 
X-ray luminosity, L x for 7 nearby Seyfert galaxies (red crosses). 
The fitted correlation W A oc I" 056 * 017 (p Spemm . dn = -0.75, sig- 
nificance 0.05), is shown as red, dashed-dotted line. When in- 
cluding NGC 253 (blue circle) - despite its questionable AGN 
nature - a correlation Wa °c ^ x °- 70±013 (p Speannan = -0.83, sig- 
nificance 0.01) is found (blue dashed line). An upper limit for 
Circinus (gray arrow) is given. The non-detection of [Ne n] from 
the nuclear point source is probably due to significant obscura- 
tion by dust in the host galaxy or the circumnuclear dust torus 
(see Fig. H]i. 



shape of continuum emission underlying the emission line. This 
uncertainty is in part due to the presence of sky lines in those 
spectra which were taken under unfavorable atmospheric con- 
ditions (e.g. high amount of precipitable water vapor). This af- 
fects the [Ne n] line the most since strong atmospheric bands are 
present between 12.4 - 12.7 yum. On the other hand, the error 
bars of the [Ar m] line are predominantly caused by the general 
weakness of the line with respect to the continuum. 

In addition to the line Wa&, Table Q] lists (absorption- 
corrected) 2 - lOkeV X-ray luminosities, Lx, for our sample 
of AGN. These are taken as a measure for the ionizing contin- 
uum luminosity since the UV/optical AGN emission of our type 
2 sources is heavily affected by dust absorption from the pu- 
tative circumnuclear dust torus. Unfortunately, 3 of our objects 
are Compton-thick to the X-ray emission (NGC 1068, Circinus, 
and ESO 428-G14). For these objects, we estimated the nuclear 
X-ray luminosity by extracting the 12yum continuum emission 
from our spectra, which will be presented in paper II, and trans- 
lating the corresponding spectral luminosity into Lx using the 
most recent version of the Lmir - Lx -correlation for AGN as 
estab lished from VISIR observations with similar spatial resolu- 
tion dHorst et al.ll2008l) . We adopt 50% of error on L x for these 
objects. 

In Figs.Q]&[2] we show the measured Wa of [Arm], [S iv], 
and [Ne n] versus L x . The trend of smaller Wa with luminosity is 
evident. As an illustration, the spectral regions around the [Arm] 
and [S iv] lines are presented in Fig.[3]for NGC 3227 (logLx = 
42.4), MCG-5-23-16 (log L x = 43.2), and NGC 1068 (log L x = 
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Table 1. Equivalent widths and X-ray luminosities of our sample of nearby AGN. 
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object 


type " 


D L " 


10/jm scale c 


W/,([Arm]) 


WMSrv]) 


t¥,([Nen]) 


L x 


ref. 






(Mpc) 


fee) 


(10"Vm) 


(lO^yum) 


(10-Vm) 


(10 43 erg/s) 




NGC 1068 


Slh 


14.4" 


21.4 


17.13 ±0.94 


38.37 ± 1.21 


26.74 ± 1.21 


5.4 ±2.7 


1 


ESO 428-G14 


S2 


22.3 


33.2 


110.57 ± 13.43 


235.16 ±8.74 


147.40 ± 2.44 


0.33 ±0.16 


1 


MCG-5-23-16 


Sli 


39.5 


58.8 


24.55 ±3.12 


48.04 ± 1.30 


111.61 ± 1.46 


1.7 ± 1.0 


2 


NGC 4507 


Slh 


54.5 


80.9 


21.41 ± 1.64 


49.14 ±3.03 


107.91 ± 16.32 


2.5 ± 0.8 


3 


Circinus 


Slh 


4.2" 


6.2 


24.66 ±2.31 


16.50 ± 4.20 


< 27 


0.59 ± 0.29 


1 


MCG-3-34-64 


Slh 


63.3 


94.2 


118.78 ±29.46 


227.63 ± 1.28 


74.35 ±2.77 


0.6 ± 0.3 


3 


NGC 3783 


S1.5 


44.7 


66.5 


23.15 ±4.00 


54.56 ±3.33 


68.99 ± 13.28 


2.0 ± 1.0 


3 


NGC 3227 


S1.5 


20.4 


30.4 


109.48 ± 18.02 


175.37 ±5.09 


275.10 ±32.69 


0.25 ±0.15 


4 


NGC 4593 


S1.0 


42.0 


62.5 


61.52 ±7.19 






0.83 ± 0.50 


5 


NGC 253 


S 


3.5 f 


5.2 






842.27 ±11.13 


0.08 ± 0.04 


1 



— Notes: " AGN types from lVeron-Ce ttv & Verq3 d2006h : * Lumino sity distance ba sed on CMB reference fr ame redshifts from N ED, other- 
wise indicated; 1 Diffraction limit for a 8.2 m VLT telescope ~0T 25; rf lTullvl dl988T) : e lFreeman etafl dl977T) : / lRekola et all j2005h 

— References: (1) L M , P -U -correlation fr omlHorst et alj(l2008h usi ng 12^m VISIR flu xes for NGC 1068 (14 Jy), Circ inus (10 Jy), ESQ 428- 
G14 (0.22 J y), and NGC 2 53 (2.1 Jy)); (2)|Balestra et aL U2004h ; (3) lHorst et ail d2008h; (4) mean of r ecent 2-10 keV dCappi et alj|2006l) and 
17-60 keV ( Sazanov et al.ir2 007) data taken, 60% error inferred due to variability; (5) Shinoza ki et alj d2006l) 



43.7). A statistical analysis of the correlation for our sample 
shows Spearman ranks pspearman = -0.75 ... - 0.9 with signif- 
icance levels in the range of 0.2 - 5 x 10~ 2 . Despite the limited 
number of objects and luminosity coverage, correlations evolv- 
ing out of the data can be confidently established. For BLBE 
and oNLBE studies, higher significance is achieved by aver- 
aging spectra of different objects within a narrow luminosity 
range. This overco mes peculiarities of individual objects (e.g. 
ICroom et al.ll2002l) . We aim for observations of a larger sample 
to make similar studies for the presented iNLBE. It is impor- 
tant to note that no correlation of Wa with AGN distance, Di, is 
present in our sample (pspearman < 0.25, significance > 0.6 for all 
lines). 

As of yet, only one other study mentions a possible 
Baldwin effect in the m id-infrared. In an AAS abstract, 
iKeremedjiev & Haol(l2006l) presented the detection of a Baldwin 
effect for the [S iv] line in AGN data obtained by the Spitzer 
satellite. They also note indications of a Baldwin effect for 
[Nen], admitting that their study is suffering from the low spa- 
tial resolution of the Spitzer data that they used. In the available 
abstract, nothing is mentioned about a slope or the scatter of the 
anti-correlation. Here, we demonstrate that the Baldwin effect of 
both the [Siv] and [Nen] line, and in addition the [Arm] line, 
is quite significant - and strong - when using high spatial reso- 
lution data, even with a small object sample. Thus, high spatial 
resolution appears to be crucial. 

3.2. Relation to M B h and fedd 

For the three type 1 AGN in our sample, black hole masses, 
Mbh, have been estimated based on reverberation mapping data. 
We adopt black hole masses for NGC 3227 (log M BH (M©) = 
7.63 + 0.31) and N GC 3783 (logM BH = 7.47 + 0.08) from 
lOnken et alj d2004l) and logM RH = 6.99 ± 0.10 for NGC 4593 
from iDennev et al" d2006l) . In addition, the black hole mass of 
the Seyfert 2 galaxy NGC 1068 could be determined by MA SER 
cloud kinematics as log Mbh ~ 7.0 dGreenhill et al.l fl996). For 
these four AGN, we analyzed the dependence of the equivalent 
widths of the [Arm] line on Mbh- A statistical test shows no ev- 
ident correlation for our limited sample of objects and the small 
coverage of black hole masses (pspearman = 0.40, significance 
0.6, for a nominal fit W / i([Arm])oc M^ ±0J ). Subject to the limi- 
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Fig. 3. [Ar m] (left) and [S iv] (right) spectral region of 
NGC 3227 (logLx = 42.4; solid blue line), MCG-5-23-16 
(logLx = 43.2; dashed green line), and NGC 1068 (logL x = 
43.7; dash-dotted red line) illustrating the decrease of Wa with 
X-ray luminosity. The spectra have been smoothed to flatten the 
continuum. 



tations, this indicates a fundamental difference of the iNLBE as 
compared t o the BLBE for wh ich such an anti-correlation has 
been found dWarner et al.ll2004l) . 

By using bolo metric correction t o the X-ray luminosities 
listed in Table [TJ( Marconi et al. 2004), it is possible to estimate 
the Eddington ratio, fedd, of four AGN with known Mbh- A sta- 
tistical test for a correlation between W^([Ar in]) and /ndd reveals 
a Spearman rank pspearman = -0.8 with a significance of 0.2 for 
the relation Wi([Arm])oc / E °j 40±0 ' 17 . Thus, our limited sample 
doesn't allow us to firmly establish a correlation, but a negative 
dependence of Wa on / E dd is indicated. 

3.3. The case of Circinus 

We excluded Circinus from all analysis of the correlations. As 
can be seen in Figs. [T]&[2] the W/s of all three mid-infrared 
lines deviate significantly from the nominal fit to the other AGN. 
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Fig. 4. Relative deviation of the observed from the correla- 
tion for all three lines plotted against the observed optical depth, 
Tgj^m, in the silicate feature. Circinus, which is an outlier in 
the Baldwin correlations of all three lines, is the object with the 
deepest silicate feature. Due to radiative transfer effects, the ac- 
tual obscuration of the central AGN by dust from the putative 
torus might be even stronger than suggested by the t^ j ^ m value. 

Although it could be a resolution effect (Circinus is the near- 
est AGN in the sample), NGC 1068 has about the same spa- 
tial resolution when scaled for the different luminosities (scal- 
ing r oc L 1 / 2 ). Moreover, we did not detect extended emission 
in the lines down to ~5% of the peak flux. This gives us some 
confidence to reject that a significant amount of line flux is over- 
resolved and lost in the noise. 

Circinus is known for strong dust extinction towards the nu- 
cleus by dust lanes in both the host and our own Galaxy. In 
Fig. |U we plot the deviation from the nominal correlation with 
X-ray luminosity (as given in Figs.Q]&[2]) of each line against the 
depth of the silicate feature. Circinus shows the strongest devia- 
tions from all correlations and, at the same time, has the deepest 
silicate feature. It is deeper than in any other type 2 AGN in our 
sample. In case of foreground extinction (host or our galaxy), we 
expect that the flux within the lines is reduced in the same way 
without a change in W,\. On the other hand, if part or all of the 
line emission is originating from inward of the torus, the extinc- 
tion in the lines could be much stronger in a more or less edge-on 
torus geometry. Given the fact that the bulk of line emission in 
all other AGN is coming from the unresolved point source, such 
a scale and inclination effect appears possible. Strong extinction 
towards the N LR based on mid-i nfrared spectra of Circinus was 
also noted bv lRoche etaT] d2006l) . 

4. Discussion and conclusions 

4.1. A comparison to UV/optical lines 

The negative correlations in the [Arm], [Siv], and [Nen] lines 
as observed with VLT/VISIR is much stronger than any known 
Baldwin effect in optical/UV broad or narrow lines. This is re- 
markable since the spatial resolution in the optical of the most 
recent studies should be similar or even better than what is 



achieved in the mid-infrared with the VLT. In UV/optical lines, 
a trend for larger correla tion slopes with incr easing ionization 
potential is observed (e.g. lDietrich et al.ll2.002l) . The potential of 
the presented lines is in the range of 40^-7 eV and, thus, com- 
parable to some of the optical/UV lines. From the optical/UV 
lines, slopes of the order of -0.1 to -0.2 would be expected. 
The measured slope of the mid-infrared lines deviates by at least 
4-cr from that expectation (except [Ne n] w/o NGC 253: 2-cr). 

In Sect 13.21 we found some negative trend of ^([Arin]) 
with /Edd, but possibly not with Mbh- If confirmed, this is in con- 
trast to the obse rved and theoretically-expected b ehavior of the 
broad lines (e.g. lWandelll999UWarner et al.l2004 . Thus, the un- 
settled correlation with Mbh might imply that the iNLBE has a 
different physical origin than the BLBE, in particular question- 
ing the relevance of the MfiH-dependence of the ionising contin- 
uum on the mid-infrared narrow-line emission. 

4.2. Implications for the origin of the iNLBE 

One suggestion for the origin of the oNLBE relates to the size 
scaling of the narrow-line region (NLR), for whi ch a luminos- 
ity d ependence ^nlr 00 L 05 has been found dBennert et al.l 
2002). As a consequence, the NLR of luminous sources may 
reach galactic size scales and, thus, lose its gas. This scenario 
is often referred to as th e "disappearing NLR" ( Croo met al.l 
120021: iNetzeretal] 120061) . For the mid-infrared lines, we find 
that the main portion of their emission is not coming from a 
spatially-extended region, e.g. as for [Om] in the optical, but 
from scales smaller than lOOpc. If the lack of line emission in 
Circinus is really an inclination-dependent obscuration effect, 
then the involved scales are significantly smaller. Since it is dif- 
ficult to imagine how the NLR could disappear on these scales, 
the iNLBE probably requires a different explanation. This ex- 
planation needs to explain, in particular, (1) the steepness of the 
correlation of the presented lines, and (2) the fact that the mid- 
infrared lines do not follow the ionisation potential-dependence 
as observed in optical lines. Possible mechanisms might involve 
a luminosity-dependent increase of density in the inner NLR 
(thus surpressing forbidden line emission). 
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